EVALUATING WHITTAKER FUNCTIONS AND MAASS FORMS 

FOR SL{3,Z) 



BORISLAV MEZHERICHER 

Abstract. We present and compare several algorithms for evaluating Jacquet's 
Whittaker functions for 5L(3,Z). The most suitable algorithm is then applied 
to the problem of evaluating a Maass form for SL{3, Z) with known eigenvalues 
and Fourier coefficients. 



1. Introduction 

Our goal is to provide tools for numerical experimentation with Maass forms for 
5'L(3, Z). The hope is that the information obtained from these experiments will be 
useful in formulating and testing hypotheses, which later can be rigorously proved. 
Even in the absence of such rigorous proofs, the data obtained by experimentation 
can lead to useful insights and intuition. 

The main computational tool is an explicit version of the Fourier expansion 
of a Maass form for SL{3,Z), given by Goldfeld in IG0IO6I . This expansion, de- 
rived in greater generality by Piatetski-Shapiro in |PS75] and independently by 
Shalika in |Sha74] . utilizes generalized Whittaker functions introduced by Jacquet 
in [Jac67] . Evaluating these functions efficiently presents a major difficulty in eval- 
uating Maass forms. Surprisingly, with the exception of the work of K. Broughan 
|Bro09j where he computes Jacquet's Whittaker functions for GL{n, R) using Stade's 
integral representation from |Sta90] , we are not aware of any computational results 
in that area. We address that problem in Section [3] where we present several algo- 
rithms for evaluating Jacquet's Whittaker function for SL{3,Z). One algorithm is 
determined to be particularly suited for evaluating a Maass form using its Fourier 
expansion. 

A possible reason for the scarcity of results related to evaluating higher-rank 
Whittaker functions is that, until recently, there were no explicitly known examples 
of generic Maass forms for SL{3,Z). This has changed last year: first examples of 
i-functions of degree 3 were discovered by C. Bian and A. Booker |Boo08] . |Bia09] 
followed by hundreds of other examples discovered by D. Farmer, S. Koutsoliotas 
and S. Lemurell |FKL08| . According to the so-called converse theorem, the Dirichlet 
coefficients of each of these L- functions should correspond to the Fourier coefficients 
of a Maass for SL{3,Z). 

With these recent developments, an algorithm for evaluating a Maass form given 
the coefficients would be very useful for further experiments. One such experiment 
is given in Section |4] where we evaluate the Maass form that corresponds to one of 
the L-functions found by Bian and provide informal evidence of its automorphy. 
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2. Preliminaries 



We need to recall some facts concerning Barnes-type integrals, automorphic 
forms for 5L(3,Z), and numerical integration. 

2.1. Barnes-type integrals. Most functions that are of interest to us can be 
represented by a Barnes-type integral, that is, an inverse Mellin transform of ratios 
of gamma functions: 



(1) m 
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1 — fem+l 
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where we adopt the notation 
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The path of integration is taken to be a vertical line from (a — ioo) to (cr + zoo), 
possibly indented so that all poles of r(l — aj + s), j = 1,2, ... ,n lie to the left of 
the path, and all poles of T{bi — s), i = 1, 2, . . . , m lie to the right of the path. For a 
description of other possible contours of integration and corresponding conditions 
on Qj's and bi's for the integral to converge, see |Luk69| (where such an integral is 
called a G-function). 

One special case is the K-Bessel function, defined below (compare to equation 
(11) in section 6.4 of [L^k69]): 

1 
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By moving the line of integration to the right, we can see that the K-Bessel function 
decays rapidly as y — > 00. Also, it satisfies 

((27ry)2 + (2^y) _ 27TyK'^ {27ry) 



(3) 



K {27Ty) 



{2^yY 

as can be seen from the definition Q by differentiating under the integral sign and 
applying the well-known recursion formula for the gamma function 

(4) r(s + i) = sr(s). 

In certain Barnes-type integral ([T|) can be evaluated by moving the line 

of integration to the left or to the right and summing the residues, using the fact 
that the gamma function has poles at non-positive integers with residues given by 

(-1)" 



Res(r(s), s — ~n) 



(n = 0,l,2,...). 



When working with the resulting series, it is convenient to use the Pochhammer 
symbol defined as 

T{x + n) 



(5) 
Note that 
(6) 
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2.2. Automorphic forms. We recall some general facts from the theory of auto- 
morphic forms on GL{n,M.); these facts can be found in |Gol06| . 
Define the generalized upper half-plane as 

l}^ ^ GL{3,R)/ {0{3,R) ■R'') . 

By Iwasawa decomposition, every z e f)'^ can be uniquely written as z = XY with 





(7) X = 1 and y = 



where Xi,yi G M and yi > 0. The group SL{3,Z) acts on f)^ by matrix multiplica- 
tion, and we are interested in functions defined on t)'^ invariant under this action. 

Let ly = (vi, V2) G C^. We introduce the following three parameters {Langlands 
parameters) 

a = -vi -2^2 + 1, 
(8) /3 = 2(^1 + 1/2-1, 

7 = -1^1 + 1^2- 

Notice that a + jS + -f ^ 0. We will also assume that ^{a) = 3fJ(/3) = 3^(7) = 0, 
that is, we assume Selberg's eigenvalue conjecture (see |Gol06j . Conjecture 12.4.4). 
We also define 

Ai = — 1 — /37 — 7a — a/3, 
A2 = -a/?7. 

Let Z? be a differential operator on t)"^ that is invariant under the action of 
5'L(3,Z). Any such operator can be written as a polynomial in Ai and A2, where 
Ai and A2 are certain differential operators given in Chapter 6 of |Gol06| or Chapter 
II of |Bum84j . 

A Maass form of type v for SL{'i,'L) is a smooth function f{z) on which is 

(a) Automorphic: f{gz) — f{z) for all g G SL{3,'Z). 

(b) Eigenfunction: Aif{z) = \if{z) for i = 1,2 and A; as in ([9|. 

(c) Square-integrable: /^^^^ ^^^1^3 |/(z)p(i*z < 00. 

(d) Cusp form: f{z) dxi dx2 dx^} = 0. 

We have the following Fourier expansion for a Maass form /(z) (see |Gol06| . 
Theorem 6.5.7): 








(10) _ A{mi,m2) 27ri[mi(c£C3+da;i)+m2SRf|^] 

^—^ ^—^ ImiTTJnl 



(c,rf) = l mi = l ms^O ' ' 



X W[ miyi\cz2 + d\, ■ 



\cz2 + d\\ 

where A(mi, 7712) G C are the Fourier coefficients, Z2 — X2 + iy2 and a, 6 G Z are de- 
fined by ad— be = 1, and M^(yi, 2/2) is a special case of Jacquet's Whittaker function. 
These higher rank Whittaker functions were introduced by Jacquet |Jac67j in much 
greater generality; for our purposes it suffices to define Jacquet's Whittaker function 
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by the following explicit integral representation initially derived by Vinogradov and 
Takhtadzhyan |VT78] and later, in more generality, by Stade |Sta90j : 



(11) W{yuy2) = m7:y,y-i{Try2y+i J K ^ {2nyiVl 



du 



u 



From the decay properties of the K-Bessel function, it is clear that Jacquet's Whit- 
taker function decays rapidly for — > oo, i = 1, 2. For a more precise description 
of its asymptotic behavior see |BH95] . 

Note that we suppress from notation the dependence of W{yi,y2) and f{z) on 
the Langlands parameters (a,/?, 7). Moreover, our definition of 14^(2/1,2/2) differs 
from other definitions in the literature by a factor dependent on the values of these 
parameters. Since (a, (3, 7) remain fixed throughout much of the discussion, we 
hope that this will not cause any confusion. Sec [Bro09 for a reconciliation of 
different definitions of Jacquet's Whittaker function as well as for the "correct" 
value of the constant factor. 

A Maass form /(z) gives rise to an L-function defined by the Dirichlet series 



(12) Lf[s) - ^ 



00 



n=l 



which has analytic continuation to the entire complex plane and satisfies the fol- 
lowing functional equation: 

s ~ a\ f s ^ (3\ f s — 7 



Conversely, under certain conditions (see |Gol06| . Theorem 7.1.3) a Dirichlet series 



such as ( 12 ) satisfying ( 13 ) gives rise to a Maass form /(z) for S'L(3, Z) with Lang- 
lands parameters (a,/?, 7) and Fourier coefficients given by the following identity 
( ]Bum8 4]. Chapter 9): 



EE 



A{m,n) _ Lf{si)Lf{s2) 
m^in^2 C{si + S2) 



The only explicitly known examples of Maass forms for SL{3,'Z) are the so- 
called lifts: given a Maass form for S'L(2,Z), a self-dual form for SL{3,'Z) can be 
constructed via what's known as Gelbart- Jacquet lift (see [G0IO6. for details). The 
generic forms, although more numerous, have been completely unknown until the 
recent numerical work by C. Bian and A. Booker |Boo08j . |Bia09| and D. Farmer, 
S. Koutsoliotas and S. Lemurell [FKLOS^ , who have produced approximations for 
the Langlands parameters and Fourier coefficients for a number of Maass forms 
by computing L-functions and using the converse theorem above. Verifying the 
automorphy of these presumptive forms is one possible application of the algorithm 
for evaluating the Whittaker functions presented in the later section. 

2.3. Numerical integration. To evaluate Jacquet's Whittaker function one can 
directly compute the defining integrals. Our method of choice for numerical inte- 
gration is essentially a version of the trapezoid rule. We now recall error estimates 
for this method from numerical analysis. 
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We apply the trapezoid rule to computing the integrals of the form fi^) dx, 
where f{x) has rapid decay as x — s- ±00. The trapezoid rule in this setting gives 

00 

^ 00 

(14) J f{x)dx^h f{kh). 

-00 fc=-oo 

It is a well-known fact (see, for example, |Fet55] . |Mor05] . or |Rub05j ) that this 
approximation is extremely accurate: the discretization error of this approximation 
is of size 0{e^'^/^). More precisely, using the Poisson summation formula, it can 
be shown that the error of the approximation above is given by X^/c^i^o / (s)' "^ti^re 
f{y) — !°"oQ 'd''" denotes the Fourier transform of f{x). If f{y) decays 

exponentially, we arrive at the error estimate given above. 

Similarly, when this method is applied to the inverse Mellin transform integral 

(15) g{y) = — j g*{s)y-^ds^- ^ g* {a + ikh)y-^-''^\ 

, , A; — — 00 

the error of this approximation is given by 

g{ye )e " . 

fc^O 

In our applications, g{y) has exponential decay as ?/ — + 00 and is bounded as 
y 0. One can see that the sum above decays exponentially if fcu < and doubly- 
exponentially if ka > 0. The resulting discretization error can be made quite small 
by choosing an appropriate value of cr. 

In our applications, the summands in the infinite sums on the right-hand side 
of the equations above decay exponentially. The sums can thus be easily trun- 
cated when several consecutive terms are smaller than the desired accuracy. The 
truncation error is then roughly of the size of the first discarded term. 

Because of the finite precision of the machine representation of real numbers, 
numerical results suffer from roundoff errors. We have done little to control this 
type of error: its effects are less pronounced when one increases working precision, 
and multiple precision arithmetic is readily available. We will point out which algo- 
rithms are most susceptible to roundoff error; improving those algorithms remains 
a research goal. 

We do not work out the error estimates in more detail than given above: in prac- 
tice, it sufhces to know how each parameter affects the error. One then chooses 
parameters experimentally by comparing results obtained with one set of param- 
eters to more accurate results obtained with a "better" set of parameters, or to 
values whose accuracy is known from theoretical considerations. 



3. Evaluating Jacquet's Whittaker function 



The main tools for evaluating Jacquet's Whittaker functions are Stade's for- 
mula (11) and the following representation in terms of a double inverse Mellin 



6 



BORISLAV MEZHERICHER 



transform: 



(16) l^(yi,y2) 



(27ri)2 



('^l)('^2) 



S2—a S2 — I3 S2— 7 si+a si+p ^1+7 
2 ' 2 ' 2 ' 2 ' 2 ' 2 
S1+S2 
2 



X (7ryi)i-«H^y2)'-^^ rfsids2. 



The equivalence of the two representations for W{yi, 1/2) can be easily established 
by applying the double Mellin transform in yi and 2/2 to both expressions. 
Note that since (a, /3, 7) = (—a, — /3, —7), we have 

(17) W{v2,Vi)^W{vi,y2). 

We now present several algorithms for evaluating Jacquet's Whittaker function. 
An implementation of all these algorithms in PARI/GP |The05j is available on the 
author's website |Mez09] . In the implementation, to avoid underflow in numerical 
computations, we work with e'^'""''' VF(?/i, 2/2)- 



3.1. Stade's formula. We use Stade's formula ( 11 ) and make a substitution u 
e" to obtain 



(yi, 2/2) =4(2^2/1 )'"*(27ry2)'+ ^ / K^{2^yiVl + ^ 

(^27r?/2\/l + e-"^ e"^"" du. 



X Ka 



Because of the exponential decay of the K-Bessel function, the integrand decays 
doubly exponentially, and the integral can be efficiently evaluated by the trapezoid 
rule as explained above. The only drawback of this method is the large number of 
evaluations of the K-Bessel function, which could be quite costly. That drawback 
aside, we found this method to be applicable to the widest range of values of yi 
and 2/2, and used it for an informal check of results produced by other integration 
methods. 



3.2. Power series expansion near (0,0). From equation (16 1, by moving the 
line of integration to the left of the poles of the gamma functions and using the 
residue theorem (see [Bum84] for details), we get 



W^(2/i,2/2)= (^2/i)'^'n^2/2)'" 



m— n— 



82 — 8x _ 82—83 
2 ' 2 

82-81 
2 



82 — Si _ ^ »53~i?i 
2 

m — n 

(-l)(™+")(7ryi)2»(7ry2)^" 
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where the outer sum runs over the six permutations of the triple (a,/?, 7). Using 
equation ([6]), we can rewrite the expression above as follows: 



§2 - h S2 - 61 S3 - di 



W{y^,y2)= (^yi)'+'n^2/2)^-'^r 
{81.52,53) 

(1 + (1 + (1 + (1 + ™! ■ 

To evaluate the coefficients of these series, one only needs six values of the gamma 
function. The rest of the coefficients can be easily evaluated recursively and, if 
desired, stored for all later calls to the Whittaker function routine. 

Note that the series converges for all values of yi and y2- However, for large 
values of yi or y2 the roundoff error becomes too large for this method to be 
applicable. To see this, note that the summands increase when, say, yi increases, 
but the resulting sum decreases exponentially. Therefore, large terms must cancel 
out, increasing the effect of roundoff errors. 

An advantage of this approach is that, given (a,/?, 7) and yi, j/2, one can easily 
determine the number of terms needed to approximate W {yi,y2)- Once working 
precision is increased so that cancellation is no longer an issue, the results of the 
computation only suffer from truncation error, and thus can be rigorously proved 
to be correct to any desired accuracy. 



3.3. Power series expansion for a small argument. In this section, inspired 
by the methods of |BSV06| . we develop an alternative power series expansion which 
provides a very good way of evaluating Jacquet's Whittaker function whenever one 
of the arguments is small. 



From equation (16 1, by moving the line of integration in si to the left and using 

(-l)"(^yi)2" 



the residue theorem, we have 

00 

W{y^,y2) =(^yi)i+"^r 



n=0 



/3 — a 7 — a 
n, n 



1 



X — / r 

27ri 



2 ' 2 ' 2 

S2 — a 



— n 



iny2r-'' ds2 



a — l3 7 ~ /3 

n, n 

2 ' 2 



(18) 



{ny,y+''Y,T 



2Tn J 

a — 7 



S2— g S2—1 
2 ' 2 ' 2 

S2-I3 



n 



n=0 



/3-7 
n, — n 



{ny2y-'' ds2 
(-l)"(7ryi)2« 



X / r 



2 — a S2~P S2—f 
2 ' 2 ' 2 

S2 —7 



(712/2)'-^^ dS2. 
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We focus on one of the three sums above; the same method apphes to the other 
two after a permutation of variables. Let 



In{y) 



1 

27ri 



2 ' 2 ' 2 



(Try) "ds. 



We have, after a change of variables S2 ^ ^2 ^ f E^nd an application of equation (|6]), 



n=0 



/3 — a 7 — a 
— — ^ n 



(19) 



1 

(3 — a 7 — a 



S2— Q S2—P S2~7 

2 ' 2 ' 2 



(7ry2)'"^^ 



(7ry2)^+^ £ 



(-l)"(7ryi)^"/„(y2) 
o(l+^)„(l+^)„-!' 



The integrals /„(y) satisfy a recurrence relation 



4+1 (y) 



1 



(20) 




In addition, 

(21) /o(y) 



1 

27ri 



(Try) ds = 4X3^ (27ry) 



where we used (|2| and the fact that a + /3 + 7 0. 

Using pTj ) , the recursion for /„ given in ( 20 1 , and the differential equation ([3| 
for the K-Bessel function, we can write 

/„(y) = (-2)-" {Pn{y)K^ (2^y) + 2nyQ^{y)K'^ (27ry)) , 
with = and polynomials Pn{y) and Qn{v) defined recursively by 

P„+i(y) = yP^{y) + ((2^y)2 + Q^(y) + a„P„(y) Po(2/) = 4 

Q„+i(y) = P„(y) + yQ:,(y) + a„g„(y) Qo(2/) = 0, 



where a„ = ^ + 2n + 2. 
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Substituting this result in (19), we get 



l+a 



r 



a 7 ■ 



E 

n=0 



i-^n^yi) 



2n 



(Tryi) 



l+a 



r 



1 - 

a 7 - 



a-/3 
2 



1 



a— 7 
2 



-^(2/2 



(7r2;2) 



2 ' 2 

Pn{y2)K^ (2712/2) + 27ryQ„(y2)if; (27ry2)) (^2/i)' 



n=0 



1 



a-13 
2 



1 



a— 7 



) 2 



with Pn{y), Qn{y) as above. Similar expressions are obtained for the other two 
sums in ( jTs] ). 

These polynomials P and Q can be computed once and stored for each sub- 
sequent evaluation of the function. Then the Whittaker function can be quickly 
and accurately evaluated with six calls to a K-Bessel routine and evaluation of 
polynomials. 

This method is very effective for small values of 2/2! if instead yi < 2/2 is small, 
one should apply equation (17 1 before using this algorithm. In fact, whenever 
either one of the arguments of (2/1,2/2) is of a reasonable size, we found this 
method to be the best way of evaluating Jacquet's Whittaker function. It avoids 
numerical integration and uses power series in only one variable. The problems 
with this method start to surface when both arguments are large: we run into the 
cancellation problems just like in Section |3.2[ For large arguments we resort to 
Stade's formula from Section |3TT] 



3.4. Numerical inverse Mellin transform. When one evaluates a Maass form 
for S'L(3,Z) using the Fourier expansion (10 1, the Whittaker function has argu- 
ments of a very specific form: one repeatedly computes W^|^2/ikl: jfp) for fixed 2/1, 
2/2 and varying z. Now let D = 2/^2/2 and define W*{D,y2) = 1^(2/1,2/2)- Then 



W\yM 



2/2 



2/2 



When z changes, D stays fixed, effectively making W* a single-variable function. 
We seek to exploit this fact in our implementation. 



From equation (16 1, we have 
1 



W{yi,y2) 



47r2 



r 



o-2+»t2~a a-z+it-z—p iJ2+it2~7 ui+iti+a cri+iti+P g-i+iti+7 ' 
2 ' 2 ' 2 ' 2 ' 2 ' 2 

cn+iti+a2+it2 
2 



X (7r2/i)i-^^-'*n'r2/2)'-"^-^*^ dh dt2 
(tt'^D) 2 (7r2/2) 2 



47r2 



(72 + ii2 - a 0-2 + «i2 - /? 0-2 + ^^2 - 7 



X F 



CTl+itl+g cTi+iti+Z^ g'i+iti+7 
2 ' 2 ' 2 

a'l+''tl+a-2+it2 
2 



(7r3D)^(7r2/2)-'*^+'^dtidt2. 
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After a change of variables, we obtain 



W{v^,y2) = W*{D,y2) 



(22) 



27r2 



r 



2 



o"l+'''2+i(*2+3ti) 



X r 



(72 + 2(t2 + ii) - a 0-2 + i(t2 + ii) - /3 0-2 + i(i2 + ^i) - 7 
2 ' 2 ' 2 

X (7r3i?)-**n'r2/2)"'*^diic/t2. 



We have written W^(yi,2/2) as a function of Z? = y\y2 and ?/2, and now exploit 
this representation in the implementation as follows. Given the parameters (a, 7), 
one pre-computes all the gamma factors needed to evaluate the integral in ( 22 1 by 
the trapezoid rule: 



(1/1,2/2) 



{■K^D) — {-Ky2) 



1 -2cr2+CT 



27r2 



JV2 



-/11/12 
2 



ikihi, 



ikihi 



Tl+CT2+''(*2+3fcifei) 

2 



X r 
X r 



f72 + i{k2h2 + kihi) — a (T2 + «(fc2^2 + ^'1^1) — P 
2 ' 2 

(72 + 1(^2/12 + fcl^l) - 7 



—ikihi 



(7r2/2) 



For some choice of parameters /ii , /i2 , ci , 0^2 , -^1 and A'2 this will result in a very 
good approximation, see ( 15 1 and the discussion that follows. 

Now, given a value of D, the inner sum is evaluated for each value of ^2 and the 
results are stored. For each subsequent evaluation of Jacquet's Whittaker function 
with the same value of D and varying y2, one only needs to re-evaluate the outer 
sum — a very fast operation. 

The main drawback of this method is that it is prone to roundoff errors, especially 
for larger values of 2/2- This can be seen quite easily: the Whittaker function is 
exponentially decaying in 7/2, but the integrand can be quite large, especially for 
larger values of \(3\ and |7|. Thus enormous cancellation must occur, leading 
to inevitable loss of precision. Still, we have found this approach to be the fastest 
of all the ones we explored. We increase working precision and use this method to 
evaluate Maass forms for SL{3,'Z) in the next section. 



4. Evaluating Maass forms for SL{3,Z) 

Suppose we are given the Langlands parameters (a,/?, 7) and the Fourier coeffi- 
cients A{mi,m2) of a Maass form f{z) for S'L(3,Z). The problem is to evaluate it 
at some point zq G i)^. 
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Since every Maass form for SL{3,Z) is even (see |Gol06j ). we can rewrite equa- 
tion ( 10 1 as 



^1 X2 /yiy2 0^ 

/(z) = / I I 1 a;i yi 
,0 1 / \ 1; 



OO OO 



= E 



(23) 



mi — 1 7712 — 1 



A{mi,m2) 

171x1712 



COS (2^1712X2) COS (2TrmxXi)W {7niyi,i7i2y2) 



Ecos (27rmi(ca;3 + dxi)) cos ( 27r— (a — 3? 



c>l 
(c,d) = l 



CZ2 + C? 



X W[ miyx\cz2 + d\, . — — 

where Z2 = 2:2 + ?i/2, a € Z is defined by ad = 1 (mod c) and (2/1,2/2) denotes 
Jacquet's Whittaker function. 

Now, given (a,(3,j) and an accuracy goal e, one can determine C — C{a^[3,^) 
such that |W^(?/i, 2/2)1 < e whenever 2/1 > C or 2/2 > C- For example, C can be 
easily determined numerically by evaluating the Whittaker function for increasing 



values of 2/1. Once C is known, one can truncate the inner sum in (23): for a fixed 
pair c, d, a term contributes whenever 



Wi2/i|cz2 + d\ < C and 



"T-2^/2 
|CZ2 + dp 



< C. 



Therefore, c and d must satisfy 



(24) 



— — < \CZ2 + d\< . 

C TO12/1 



It follows that for each miyi, m22/2 and Z2, the infinite sum over all relatively 



prime integers c and d with c > 1 can be truncated whenever the condition (24) 



is not satisfied. Because of the exponential decay of Jacquet's Whittaker function, 
truncation of the outer sums over mi and m2 does not pose a problem. 

To evaluate the inner sum for fixed values of mi2/i, TO2J/2; and Z2, we use the 

D — (TOi2/i)^m22/2 is fixed, and we compute the value 
for all c, d that satisfy (24). 



algorithm from Section 
of W I mi2/i|cz2 + d\ 



3.4 



|cZ2+(i|^ 



As mentioned in Sec- 



tion when D is fixed, each Whittaker function evaluation takes very little time. 
The only time-consuming operation is changing the value of D, but that happens 
much less frequently during the computation. 

We apply this algorithm to evaluate several 5'L(3, Z)-Maass forms f{z) at zq = 
/s, the 3x3 identity matrix. The Maass forms are lifts of Maass forms for SL{2, Z), 
which have been computed to high precision in |BSV06] . The Langlands parameters 
of lifted forms are (— 2ir, 2ir, 0), where r is the type of the SL{2, Z)-Maass form as 
in |BSV06] . The time required for each computation and the number of the Fourier 
coefficients needed to evaluate /(zq) with error not exceeding 10~^^ are summarized 
in the table below. 

The number of coefficients is essentially the maximum value of m2 for which we 
get a non-negligible contribution from the innermost sum in ( 23 1 . As should be 
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BORISLAV MEZHERICHER 



r 


Cocffs 


Time (s) 


9.533695 


120 


94 


13.779751 


203 


204 


17.738563 


379 


440 


22.194674 


590 


972 


26.056918 


861 




35.431665 


1806 





Table 1 . Time and number of the Fourier coefficients needed to 
evaluate a lifted Maass form on a Pentium D 2.80GHz. 



expected from (24 1, the maximum value of mi resulting in non-negligible contribu- 
tion is much smaller, not exceeding 10 for the first cases. It would be interesting to 
see if there is another way to evaluate a Maass form that requires fewer coefficients 
for the same level of precision, perhaps by making the expansion more symmetric 
in TOi and m2. Given that Farmer et al |FKL0 8 require only about 20 coefficients 
to locate a candidate Maass form, it is somewhat surprising that evaluating such a 
form would require an order of magnitude more coefficients. 

In another application of this algorithm, we evaluate a generic Maass form with 
Langlands parameters (-14.141638,-2.380388,16.522027) discovered by Ce Bian 
|Bia09] . Let 



= -1 and 52 





be two elements of SL{i, Z). We pick several points zq = {{xi,X2^ xz), (2/1,2/2)) S f)"^ 
and evaluate /(^o) and f{'jzo) where 7 is a word in S^'s. Some results are given in 
the table below. 





f{zo) = f{xi,X2,X3,yi,y2) 


((0,0,0), (0.9, 0.9)) 
5i525i((0,0,0),(0.9,0.9)) 


-79.779900 - 0.000044759125i 
-79.780569 - 0.000025792266i 


((0,0,0), (0.9, 10/9)) 
525i((0, 0,0), (0.9, 10/9)) 


-168.098614 -H7.88405197i 
-168.099135 -l-7.88392995i 


((0.2, 0.4, -0.1), (0.95, 0.9)) 
S'25'i((0.2, 0.4, -0.1), (0.95, 0.9)) 


-16.1067874 -M2.9212510i 
-16.1066031 + 12.9212690i 


((0.726, 0.325, 0.983), (0.541, 0.578)) 
5*2 5*1 ((0.726, 0.325, 0.983), (0.541, 0.578)) 


11.3789573 + 0.0986941489i 
11.3788905 + 0.0986972163i 


Table 2. Values of a Maass form f{z) with (a,/3,7) = 
(-14.141638,-2.380388,16.522027) discovered by Bian and 
Booker 



The agreement between the values of / at two 52^(3, Z)-equivalent points of i)^ 
is quite good. Given the precision of the Langlands parameters and the Fourier 
coefficients, the errors of this size are to be expected. We thus have further evi- 
dence of automorphy of Bian's first example. A rigorous numerical verification of 
automorphy along the lines of |BSV06] remains a research goal. 
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